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ABSTRACT 

SN 2010jl was an extremely bright, Type Iln SNe which showed a significant IR excess no later 
than 90 days after explosion. We have obtained Spitzer 3.6 and 4.5 /im and JHK observations 
of SN 2010jl ^90 days post explosion. Little to no reddening in the host galaxy indicated that 
the circumstellar material lost from the progenitor must lie in a torus inclined out of the plane 
of the sky. The likely cause of the high mid-IR flux is the reprocessing of the initial flash of 
the SN by pre-existing circumstellar dust. Using a 3D Monte Carlo Radiative Transfer code, we 
have estimated that between 0.03-0.35 Mq of dust exists in a circumstellar torus around the SN 
located 6 x 10^^ cm away from the SN and inclined between 60-80° to the plane of the sky. On 
day 90, we are only seeing the illumination of approximately 5% of this torus, and expect to see 
an elevated IR flux from this material up until day ^ 450. It is likely this dust was created in 
an LBV- like mass loss event of more than 3 Mq, which is large but consistent with other LBV 
progenitors such as rj Carinae. 

Subject headings: supernovae: individual (SN 2010jl) — supernovae: general — circumstellar 
matter — dust, extinction 



1. Introduction 

SN 2010jl was the brightest core collapse supernova (CCSN) of 2010, wit h a peak unfiltered bright ness 
of 12.9 mag, and a corresponding peak absolute magnitude of r oughly -20 jNewton fc PuckettI [im^ . A 



spectrum taken on 2010 November 5 was classified as a Type Iln (jBenetti et al.ll2010f ). a particular class of 
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Type II SNe which show narrow (~100 km s ^) and often intermediate (^1000 km s ^) lines in H and He 
along with normal broad (10000 km s~^) lines due to the expansion of the ejecta (!Schlegellll990|). The narrow 
emission lines in their spectra are attributed to the ionization of the pre-existing circumstellar material (CSM) 
which has been excited by the initial flash of the supernova. Intermediate lines can arise at early times due to 
the high optical depth in the CSM, whi ch can cause multip le scatterings from thermal electrons which results 
in the broadening of the narrow lines ( Smith et al.ll2010r) . It is believed the progenitors of Type Iln SNe 



which constitute 9-10% of the total population of CCSNe ( Smith et al 



2011 



Li et al 



2011 



Smartt et al 



20091), lose 10^4 to 10"^ Mq yr"! (|Fox et al.ll201ll : iKiewe et alJboiol ) orders of magnitude more material 
than normal Type II progenitors. Pre-discovery Hubble Space Telescope (HST) images indicate that the 
progenitor of SN 2010jl w as likely a massive s tar, >30 Mq, which could have been a luminous blue variable 
(LBV) in outburst phase |Smith et al.lboilbl ). 



SN 2010jl is only the second CCSN progenitor with a mass greater than 17 Mq detected, the other being 
the LBV progenitor of SN 2005gl (jCal-Yam fc Leonardlliooi ICal-Yam et allboOTh . This is important for the 
study of dust production around CCSNe s ince the SNe from less massive progenitors (8-17 Mq) only seein 
to be producing 10"^ - 10""^ Mq of du st |Elmhamdi et al.lbooj [Sugerman etaPbood iMeikle et al 112007 : 
Kotak et al. 2009t Andrews et al. 2010 . for example). This is far less than the ~ 0.1-1 Mq of dust estimated 
per CCSNe to explain the large amounts of dust seen in galaxies of z > 6. This has led to speculation that 
more massive progenitors existe d among the Population HI stars of high-z galaxies, which may produce more 
dust IjCherchneff fc Dwekll2010l . for example). To date though, there is indirect evide nce that some nearb y 
dust producing CCSNe have progenitors massive than 17 Mq, for example SN 2007it (jAndrews et al.ll2011l ). 
and are also only producing dust masses of 10^'^ Mq. This may indicate that the amount of dust produced 
is independent of progenitor mass. The study of dust formation in SN 2010jl, which likely has a massive 
progenitor, will provide an opportunity to explore this relationship more fully. 

Due to the nature of Type Iln SNe, we expect pre-existing CSM gas and dust to exist around the SN. 
Therefore to distinguish newly formed dust from dust that was already present, it is imperative to take 
early-time observations to measure the contribution from pre-existing dust to serve as a baseline for later 
times when conditions may be favorable for new dust condensation. Studying the pre-existing dust around 
CCSNe can also be essential to understanding the progenitor properties. Continued observations of the 
interaction between the SN and the surrounding CSM material over the first few years can reveal valuable 
information about densities , velocities, a,iid oy er all mass-loss history of the progenitor as well as progenitor 
mass itself ( Fox et al. 2011 : Smith et al. 20091 and references therein). This paper on SN 2010jl presents a 
spectral energy distribution (SED) using optical, near- and mid- infrared photometry of SN 2010jl obtained 
~ 90 days after explosion and estimates of the mass, size, and age of the circumstellar dust. 



Observations 



Although the discovery date of SN 2010jl was 2010 November 3.5 (Newton fc Pucketl 2010[ ). compar- 
i son of the optica l spectra with other similar SNe ( Yamanaka et al. 2010) as well as pre-discovery images 
dStoU et al-lboiol) indicate that the explosion date was likel y early October. An ASAS image taken on 2010 
October 9.6 shows the SN with a V mag nitude of 13.79 (jstoh et al.llioioh . This is consistent with the 
GELATO mod eling of SN 2010jl to SN 2006tf, which indicated that it was about a month past explosion 
(lYamanaka et al. 2010), and has led us to adopt 2010 October 10 (JD 2455480) as the date of explosion 
throughout this paper. It is located in the interacting g alaxy UGC 5189A, for which throughout this paper 
we adopt the distance of 48.9 Mpc (ISmith et al.ll2011bl ). We have also assumed a Galactic foreground red- 
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dening of E(B-V) = 0.027 from Schlegel et al. ( 19981) . and ar e assu ming zero reddening in th e host galaxy 
based on the previous estimates for SN 2010jl bv lPatat et al.l (|201l[ ) and lSmith et al.l (|2011bl ). 



A UVBRI photometric sequence of the UGC 5189A field (listed in Table 1 and shown in Figure [TJ was 
taken on the night of 2010 December 5 (JD 2 455504) at the KPNO 2.1m telescope. These tertiary standards 
were calculated using the same methods as [Andrews et al.l (|2010l 120111 ). The processing of all images to 
remove instrumental signatures was accomplished using the standard techniques of subtracting a median 
filtered bias frame and dividing by master twilight flat field images for each filter. Also, the fringing effects 
in the I-band exposures were removed by scaling and subtracting a master fringe frame from the program 
images. In order to transform the instrumental magnitudes to the standar d syste m, between 30 and 35 
different stars (depending on filter) were observed from the lists of iLandoltl (j2009[ ) in addition to the SN 
2010jl field. These stars were observed both near zenith and at a high airmass (~1.75) and were selected 
due to their color ranges. Thus, they provided a viable sample of stars from which to derive accurate color 
and extinction coeffecients as well as magnitude zero-points for each filter. The resulting RMS scatter of the 
residuals from the transformations revealed that the photometry for selected stars in the SN 20010jl field 
derived from this night's observations is accurate at the 1-2% level depending on filter . BVRI photome try 
was obtained either from observations from the KPNO 1.2m (Table 2), the literature ( StoU et al. 2010l) or 
via the AAVSO data server and is plotted in Figure 2. The data shown in the plot have not been corrected 
for foreground extinction. For the SED we used an observation from 2011 January 3 (JD 2455563), day 
88, submitted to the AAVSO by contributer Etienne Morelle (Observations from the AAVSO International 
Database, 2010, private communication). 

JHKs observations of SN 2010jl w ere taken on 2011 J anuary 21 (JD 245583), or day 108, with the 
WHIRC camera on the WIYN telescope (jMeixner et al ]l2008l) . The detector offers a pixel scale of 0.098"pix ^ , 
and a FOV of 3'x 3'. All images were taken in dither mode and aligned and stacked prior to fiat fielding 
and bias correction. Eleven frames of 120s were taken in the J band, 12 frames of 120s in the H band, and 
13 frames of 90s in the K band were taken of the SN 2010jl field. The transformation from instrumental 
to JHKs magnitudes was accomplished using 2MASS stars located in the frame for comparison photometry. 
Uncertainties were calculated by adding in quadrature the uncertainties associated with the 2MASS com- 
parison stars and the uncertainties in the PSF photometry. The photometry is listed in Table 3 and shown 
in Figure 3. 

The Spitzer IRAC (3.6 and 4.5 fim) images taken on 2011 January 5 (JD 2455565), corresponding to 
day 90, were mosaicked and resampled using standard MOPEX procedures to improve photometric quality. 
Pre-explosion IRAC images of UGC 5189A were also available in the Spitzer archive from 27 December 
2007 (Program 40301, PI Fazio) which were used to subtract from the SN 2010jl images to get accurate 
photometry. PSF photometry was performed using the position specific PRE images. Table 3 contains the 
measured mid-IR fluxes obtained from Spitzer. Figure 3 presents the SED of SN 2010jl on ^ day 90 (88 for 
BVRI and 108 for JHK). Statistical uncertainties presented in the plot repres ent la errors. T he data shown 
in Figure 3 have been corrected for foreground extinction of E(B-V) — 0.027 (jSchlegel et al.lll998i) . 



3. Discussion 

The SED of SN 2010jl, shown in Figure 3, shows an IR excess on day 90 which indicates the presence 
of warm dust around the SN. The IR excess could be explained by two possible scenarios, one being that 
pre-existing CSM dust existing beyond the evaporation radius is reprocessing the initial flash of the SN into 
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the infrared. Another possibility, although less likely, is that there is already dust forming in a cool dense 
shell (CDS) that has formed bet ween the forward and revers e shocks created b etween the ej ect a and the 



20081 ) and SN 2005ip ([Smith et al. 



pre-e xisting CSM gas interaction (Chevalier fc Franssonlll994t iPozzo et alj 120041 ) . SN 2006jc (jSmith et al 



200' 



i 



Fox et al 



200' 



^ formed dust at day 50 and 75 respectively, due to 



interaction with the CSM. This newly formed dust would also create an IR excess. Without optical spectra 
taken at the same epoch to look for the emergence of asymmetries, we cannot clearly distinguish between 
pre-existing and newly formed dust. 

The progenitor star of a CCSN is expected to have lost massive amounts of gas and dust prior to 
explosion, which will exist in the CSM. When the star becomes a supernova, the initial UV flash from the 
explosion wil l evaporate the dus t grai ns a certain dista nce from the SN, dependent on the initial luminosity. 



According to lDwekl (|l983l . 119851 ) . and lFox et al 

clear a cavity with a radius, levap = 6 x 10^^ cm if the CSM is carbon rich, or Tevap 



(I200i), a SN with an initial luminosity of 1 x 10^" L© will 

3 X 10" cm if the 



CSM is oxygen rich. Assuming the maximum luminosity of SN 2010jl was 3 x lO"'^ erg s~^ (Mi,=-20), this is 
roughly 10^'^ Lq, we can expect a cavity of cleared material around the SN to have a radius of approximately 
the same size. 

Using black body fits (shown in Figure 3), we made initial estimates of Ri„ (the inner radius of the 
shell), luminosity and temperature to constrain the model. These yielded a dust temperature (T^) of roughly 
750 K, ejecta temperatures (Tgj) of 7500 K , and an Ri„ = 1.3 x 10^^ cm for this epoch. The blackbody 
temperature and radius of the IR component of SN 2010jl suggests pre-existing dust, since new dust formed 
in the CDS prior to day 90 would be expected to be much hotter, ~1600 K, as was the case in SN 2006jc 
and SN 2005ip and much closer to the expansion center, since the ejecta o n day 90 would hay e traveled a 
maximum of 1.1 x 10^^ cm assuming a constant expansion of 14 000 km s~^ (jBenetti et al.ll2010l ). Any CDS 
would have to arise interior to this radius. Therefore, we believe that the IR emission most likely arises from 
pre-existing dust grains surrounding SN 2010jl heated by the initial flash. 

As mentioned above, pre-explosion IRAC images do exist for UGC 5189A, but we are unable to detect 
the CSM shell surrounding the progenitor. This is not surprising if we assume either a red supergiant (RSG, 
L = 5 X 10"* Lq, T = 3500K) or an LBV (L = 1 x 10^ L©, T = 10000 K) as a progenitor star, 48.9 Mpc 
away with dust at a distance of 6 x 10^^ cm, the total fiux contribution from both the star and the dust in 
the pre-explosion image in 3.6 and 4.6 fim is less than 0.5 /xJy. 

In order to estimate the pre-existing dust mass in SN 2010jl, we used our radiative transfer code 
MOCASSIN, a three-dimensional Monte Carlo radiative-transfer code that is capa ble of modeling non- 
spher i cal ge ometries, clumpy density distributions, and non-central energy sources (jErcolano et al.l 12003 . 



2005 



200 



As an initial fit, we used a "smooth" model, in which the dust is uniformly distributed 
throughout a spherical shell according to a r~^ density profile, and constrained within a shell with size R,;n 
and Rout- We used a standard MRN grain size d istribution of a~^-^ b etween 0.005 and 0. 05 u m (IMathis et al 



1977), as was used in previous modeling done by Meikle et al. (2007), Kotak et al. ( 20091) . and Andrews et al 



(j201l[ ). Due to the early epoch of the observations, we assume the SN is a point source at the center of a 
spherical shell. 

We us ed version 2.02 . 67 of MOCASSIN, which includes the capacity to account for light travel time 
effects (see IWesson et all (l20ld ) for details), to calculate the emission from only that portion of the dust 
illuminated after 90 days (Figure SJ. When light travel times are taken into account, after 90 days only 
about 5% of a spherical shell with Ri„ = 6 x 10^^ cm and Rout = 1-4 x 10^^ cm would be illuminated as 
seen from Earth (see Figure 4). Applying this constraint, spherical dust geometries could provide a good 
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match to the observed day 90 SED, with a total dust mass of ^0.6 M(7i. but g ave a Hne-of-siKht optica l 
depth of ^^6. This was problematic since previous papers from lPatat et al.l (|201ir ) and (jSmith et al.ll2011bl ) 
estimated a Galactic foreground reddening of E(B-V) = 0.027 for SN 2010jl and no reddening f rom the host 



galaxy . SN 2010.il is a goo d match to SN 2006tf, which also had little host galaxy reddening (jSmith et al 



2011bl ). iPatat et al.1 (j201ll) , using both the continuum of the optical spectra and equivalent widths of the Na 
I D lines, also found no host galaxy reddening. Therefore, a geometry other than a spherical shell, such as 
an inclined torus or bi-polar lobes, needs to be considered in order to allow low optical depth along the line 
of sight, but with sufficient dust mass to account for a CSM envelope with the high IR excess of SN 2010jl. 

For the torus models we used a uniform dust d ensity with a distanc e of 1.2 light years (6 x 10^*" cm) 
and a tube radius of 0.5 light years (4.7 x 10^^ cm) (jErcolano et al.ll2007n . A torus with zero inclination (in 
plane of the sky) would not have been illuminated by the supernova by day 90. An inclination of at least 45° 
is required for IR emission to be present on day 90. Figure [3l shows the fits for the 60° inclination. We used 
both 100% amorphous carbon (AC) compositions and 100% silicate (Si) compositions, since the lack of data 
beyond 4.5 /.tm, makes it impossible to differentiate between silicate and AC grains. We find an inner radius 
of 6 X 10^^ cm also provided the best fit to the data, which is actually larger (particularly for a carbon rich 
CSM) than the predicted evaporation radius discussed above. This suggests that dust was absent interior to 
Tetjap before the explosion, and that the pre-existing dust was actually formed by a progenitor eruption as 
early a s 300 years prior ( if we assume an LBV eruption of ^ 600 km s~^, as was seen in the bipolar lobes of 
rj Car (jSmith et al.l 120031 )) or as late as 20,000 years if we assume an RSG progenitor with a constant wind 
of 10 km s~^. 

The 45° inclined torus, although having a zero optical depth underestimates the IR flux. The 80° torus 
does a fairly good job fitting the data, but has an optical depth of 1.53 and 2.4 for the AC and Si models 
respectively, much too high for the reddening suggested by the early time optical spectra. The 60° torus 
includes both a zero optical depth and an acceptable fit to the data, making it the most likely inclination 
of the surrounding CSM torus. With only the two mid-IR points the exact inclination is not very well 
constrained, but an inclination of 60 - 80° to the plane of the sky surrounding SN 2010jl is consistent with 
the data. This model (see input parameters listed in Tables 4 and 5), implies that between 0.03 and 0.35 
M0 of pre-existing dust surrounds SN 2010jl. This would indicate total masses of CSM between 3- 35 Mq, 
assuming a gas to dust ratio of 100 and the possible inclinations and dust compositions. Although this is a 
large amount, measurements of d ust around Eta Carinae, an extensively studied LBV in our galaxy, yield 
dust values between 0.4 - 0.7 M©, ( Gomez fNee Morgan) et al. 2006 : Gomez et al. 2010l) making it plausible 
that previous mass loss events could have ejected this much dust around SN 2010jl. Like the geometry 
of rj Carinae, we cannot rule out the dust being contained in bipolar lobes and not in a torus. However, 
observations at later epochs will provide stronger constraints on the geometry of the dusty CSM which we 
have shown to exist around SN 2010jl. 



4. Summary 



This paper presents the results of SED modeling of SN 2010jl ^ 90 days post-explosion. Using optical, 
near- and mid-IR observations of SN 2010jl from day 90 we have estimated a pre-existing dust mass between 
Md = 0.03 - 0.35 Mq, with the higher dust masses corresponding to an oxygen-rich CSM. This dust is 
likely located in a torus inclined between 60 - 80° so me 6 x 10^^ cm away fro m the SN. The inclined torus 
is required by the conclusion of iPatat et al.l (|201lh and (jSmith et al.ll2011bl ) that there is very little local 
reddening around SN 2010jl, while at the same time allowing the existence of substantial amounts of dust 
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needed to create the observed IR excess. The lack of dust between the estimated evaporation radius and the 
inner dust radius suggests that this torus of pre-existing dust was created in a mass loss episode likely 300 - 
2000 years prior, assuming an LBV progenitor with expansion velocities between 100-600 km s"^. As many 
LBVs can have bipolar or toroidal nebulae and m assive dust shells, it seems very likely that the progenitor 
was an LBV. This agrees with the conclusions of ISmith et al.l (|2011bl ). who suggests the progenitor was a 
massive LBV star. If the inner edge of the CSM dust is 6 x 10"'^'' cm away we estimate that the contribution 
to the elevated IR flux from pre-existing dust will last at least for 1.2 years, and that it in fact will increase 
as more and more of the torus is illuminated. Continued monitoring of this object is crucial in order to 
detect the formation of new dust, and to measure how much dust is created. 
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Table 1: Tertiary BVRI Standards for UGC 5189 



Star U B V R I 

20.553 ± 0.1351 19.611 ± 0.0453 18.486 ± 0.0180 17.765 ± 0.0148 17.152 ± 0.0366 

B 19.087 ± 0.0429 18.468 ± 0.0176 17.513 ± 0.0123 16.963 ± 0.0084 16.488 ± 0.0200 

C 14.929 ± 0.0138 14.677 ± 0.0036 13.944 ± 0.0049 13.521 ± 0.0048 13.131 ± 0.0038 

D 19.807 ± 0.3422 19.470 ± 0.0508 18.164 ± 0.0138 17.313 ± 0.0167 16.606 ± 0.0153 

E 19.967 ± 0.0636 18.751 ± 0.0537 17.189 ± 0.0065 16.062 ± 0.0111 14.685 ± 0.0052 

F 19.888 ± 0.0995 19.456 ± 0.0448 18.212 ± 0.0201 17.511 ± 0.0154 16.823 ± 0.0450 

G 16.986 ± 0.0157 16.323 ± 0.0060 15.405 ± 0.0027 14.877 ± 0.0059 14.427 ± 0.0055 

H 18.965 ± 0.0834 19.022 ± 0.0597 18.310 ± 0.0162 17.848± 0.0220 17.504 ± 0.0274 

I 20.439 ± 0.3093 19.566 ± 0.1121 18.433 ± 0.0174 17.685 ± 0.0211 17.053 ± 0.0334 



Table 2: BVRI observations of SN 2010jl 



JD 


Age 


B 


V 


R 


I 


Telescope /Instrument 


2455535 


55 


14.32 


13.94 


13.51 


13.32 


KPN0/1.2m 






± 0.01 


± 0.01 


± 0.01 


± 0.01 




2455629 


149 


15.16 


14.69 


14.07 


14.08 


KPN0/1.2m 






± 0.01 


± 0.02 


± 0.01 


± 0.01 




2455656 


176 


15.20 


14.76 


14.09 


14.16 


KPN0/1.2m 






± 0.03 


± 0.01 


± 0.01 


± 0.01 





Table 3: IR observations of SN 2010jl 



JD Age J H Kg 3.6 /xm 4.5 /xm Telescope/Instrument 

2455565 90 4.04 mJy 4.52 mJy Spitzer/IRAC 

± 0.14 ± 0.18 

2455583 108 13.41 13.09 13.75 WIYN/WHIRC 
± 0.2 ± 0.2 ± 0.2 



Table 4: Monte Carlo Radiative Transfer Torus Models for 100% Amorphous Carbon 



Inclination 


Tei (K) 


Rin (cm) 


Rout (cm) 


Ltof. (Lq) 




Md (Mq) 


45° 


7500 


6el7 


1.4cl8 


5.5e9 


0.0 


0.12 


60° 


7500 


6cl7 


1.4cl8 


5.5c9 


0.0 


0.05 


80° 


7500 


6el7 


1.4el8 


5.5e9 


1.5 


0.03 
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Table 5: Monte Carlo Rmliative TriULsIer Torus AIcxlcLs i'or 100% Silic.aU;^ 



Inclination 


Tei (K) 


Rin (cm) 


Ro„t(cm) 


Ltof. (L©) 




Md (Mq) 


45° 


7500 


6el7 


1.4cl8 


2.6el0 


0.0 


0.7 


60° 


7500 


6el7 


1.4el8 


2.6el0 


0.0 


0.35 


80° 


7500 


6el7 


1.4el8 


2.6el0 


2.4 


0.27 




Fig. 1. — V-band finder chart for tertiary standards for the field around UGC 5189 on 5 April 2011. 
Magnitudes are listed in Table 1. 
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Fig. 2. — Optical lightcurve of SN 2010jl. Data are from AAVSO, StoU et. al 2010, and new data presented 
in Table [2] in this paper. 




Fig. 3. — SED of SN 2010jl on ~ day 90. IR photometry is listed in Table [3l and optical photometry was 
obtained from the AAVSO, specifically contributor Etienne Morelle. Dashed dotted lines are the blackbody 
fits to the data. The solid line indicates the 100% amorphous carbon inclined at 60° MOCASSIN torus 
model fit, the dotted line is the same except for the composition being 100% silicates. As the figure shows, 
without observations longward of 4.6 /im, we cannot accurately constrain the dust composition. 
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Fig. 4. — Schematic of the fraction of the CSM ejecta which wih have been ht when hght travel times are 
taken into consideration. For the example above, at 90 days, only about 5% would have been illuminated. 



